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We explore the combination of the extended dynamical mean field theory (EDMFT) with the 
GW approximation (GWA); the former sums the local contributions to the self-energies to infinite 
order in closed form and the latter handles the non-local ones to lowest order. We investigate the 
different levels of self-consistency that can be implemented within this method by comparing to the 
exact QMC solution of a finite-size model Hamiltonian. We find that using the EDMFT solution 
for the local self-energies as input to the GWA for the non-local self-energies gives the best result. 
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Introduction - The GW approximation (GWA) [l| is 
one of the most successful methods to describe elec- 
tronic structure of weakly correlated materials. It in- 
cludes the lowest order perturbative corrections in the 
screened Goulomb interaction to the electron self-energy. 
In real space-time, the electron self-energy is given by 
the product "GW" , where "G" represents the electron 
Green's function and "W" the RPA screened Coulomb 
interaction. The GWA has been successfully applied to 
the calculations of quasi-particle spectra of semiconduc- 
tors and insulators. (For recent reviews, see Ref.Q) 
It describes well the experimentally observed energy gaps 
in semiconductors 0, II, IE S 013 • 

The GWA self-energies were first calculated as a one- 
shot perturbation by using the unperturbed LDA Green's 
functions. However, as we will show in this paper, 
the one-shot GWA breaks down as soon as the correla- 
tions are moderately strong, as expected for a non-self- 
consistent leading order approximation. The problem 
should be solved by including higher order corrections 
to extend the scope of the method. This is the main mo- 
tivation of our work. The Baym-Kadanoff formulation 
of the GWA, which makes the method automatically 
conserving, requires to evaluate the self-energy, the polar- 
ization bubble, and the Green's function self-consistently 
0. While this has been successfully implemented 0| 
and shows improvement on the values of the total en- 
ergy, the self-consistent GWA brings up a debate over 
whether the fully self-consistency improves the spectra. 
II S m E 111 El- This is due to the fact that 
the vertex correction, which is omitted in GWA, shows 
a tendency to cancel the self-energy insertion when the 
full Green's function is used in calculating the exchange 
self-energy 0, 0, 0| . As a result the one-shot GWA or 
partial-self-consistent ones (in which, e.g., one fixes the 
"W" obtained from the LDA input while sol ving "G" self- 
consistently) are favored sometimes [ill IT^ Il3j |. This is- 
sue entails the usage of higher order self-energy diagrams 
and provides us a second motivation for improving sys- 
tematically the approximation beyond the leading order. 

In this paper, we extend the GWA by using the ex- 
tended dynamical mean field theory (EDMFT) [H HI 
mi^ which treats all the local polarization, self-energy 



and vertex corrections in closed form, following the ideas 
proposed in R,efs.[2(Tl I21I l2^ . According to these ap- 
proaches, the local polarization and self-energy should 
be solved by EDMFT non-perturbatively while the non- 
local ones via GWA. The cancellation of the self-energy 
insertion and the vertex correction is carried out lo- 
cally to infinite order within EDMFT. However, the 
same ambiguity of implementation (with or without self- 
consistency) remains in the non-local GWA part. This 
question can not be addressed without benchmarking the 
different schemes. 

Main Results.- We compared and contrasted the differ- 
ent implementations of combining EDMFT with GWA. 
Our numerical calculation suggests that the schemes us- 
ing the EDMFT solution of the local self-energies as in- 
put to the GWA for the non-local self-energies give the 
best result. In other words, it is favorable not to allow 
the feedback of the non-local GWA on the local EDMFT 
self-energies due to different nature of the two methods, 
GWA being perturbative and EDMFT non-perturbative. 
On the other hand, when using the EDMFT result as in- 
put, it makes little difference whether one performs GWA 
with a single shot or partial self-consistency (it is par- 
tial because the local polarization and self-ener gy f rom 
EDMFT are fixed). This supports the picture |2ll |2^ 
that, in a correlated phase far away from a phase transi- 
tion, the temporal correlations reflected through the local 
polarization and self-energy are dominant and should be 
treated non-perturbatively, while the spatial correlations 
are weaker and can be handled perturbatively. 

Lattice Model- In comparing different GW ap- 
proaches, many of which apply to realistic materials, one 
always encounters the problem that there are as many 
important differences in the implementations as in the 
GWA methodologies themselves. These include the im- 
plementations of GWA in imaginary time or frequency 
space, the choice and construction of the local basis set, 
and the further approximations like the plasmon-pole ap- 
proximation 3]. In order to compare the many-body 
schemes without additional complications, it is desired 
to employ simple model systems on the lattice 11, 23]. 
Following the same strategy, we study the following gen- 
eralized Hubbard model: 
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FIG. 1: The potential <E> of the Baym-Kadanoff functional 
[mi for (EDMFT+GW)sc. (EDMFT)sc is obtained by 
restricting the exchange diagram (the second on the r.h.s.) 
to be local in space, i = j. The first line by itself repre- 
sents the standard (GW) sc scheme. In this formulation, the 
boson Green's function D describes the screened interaction. 
The boson self-energy plays a similar role as the electron-hole 
bubble in GWA. 
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FIG. 2: In our model system, the one-shot GWA, using the 
Hartree-Fock result as input, breaks down at (7 ~ 0.825. This 
is an instability against the formation of charge density wave 
at (tTjTv) [2^. The other results are plotted as references. 
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where cia {c\^) annihilates (creates) an electron of spin 
a (=t, i) at the lattice site i. hi is the electron density. 
We only consider the paramagnetic phase. 

Approximation Schemes - The self-consistent schemes 
are summerized in FigQ] [we use {■ ■ ■)sc for a 
self-consistent (SC) loop and (• • ■)psc a partially 
SC loop.] The full functional in FigQ] can be sep- 
arated into local and non-local parts, which are 
functionals of the local and non-local Green's func- 
tions, respectively: $ = ^ EDMFt[G Local, D Local] + 
^Non~Local Gw[GNon-Local, DNon-Local]- The lo- 

cal (non-local) self-energies, which arc obtained 
from the functional derivatives of ^edmft 
N on- Local Gw), are functious of the local (non- 
local) parts of the Green's functions only: ^edmft = 

^EDMFT{GLocal, DLocal) and ^ Nan- Local GW — 
^Non-Local Gw{GNon-Local, DNon-Local) for the elec- 
tron, He DM FT — ^ED M FT {G Local, D Local) and 
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FIG. 3: The energies, the quasi-particle residue (Z), and the 
non-local contribution to the self-energy {Z\) as defined in 
the text. The same symbol scheme defined in (a) applies to 
aU the diagrams. Z is the same for (EDMFT)sc+(GW)psc, 
(EDMFT)sc+GW, and (EDMFT)sc, since the same local 
self-energy is used. For Zi, which reflects the spatial extension 
of the self-energy, we are comparing numbers at least one 
order smaller than those for other quantities. 



^Non-LocalGW — ^Non-Local GW {G Non- Local) for the 

boson. 

The approximation schemes we tested in combining 
EDMFT with GW are: (i) The fully self-consistent 
(EDMFT+GW)5c solves the full Dyson equations for 

G — G Local +GNo7i-Local and D = Dloc 
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FIG. 4: The local electron Green's function and self-energy 
at [/ = 1.5. The same symbol scheme defined in (a) applies 
to (b). The lines connecting the points are guides to the eye. 
The fermion Matsubara frequency p„ = {n//3){2n + 1). 
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(EDMFT)sc+GW uses 



The non-self-consistcnt 
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EDMFT, Non-Local — [<J"o ^ ■^EDMFTlNon-Local 
EDM FT, N on-Local — [Dq — ^EDMFT]iq'on- Local 

obtains in one shot the estimation for TiNon-Locai aw — 

^Non-Local Gw{GeDMFT .N Oil — Local ■) D EDM FT, Non-Local) 
and ^Non-Local GW — 

^Non-Local GW {G EDMFT, Non-Local ) ■ (iii) 

(EDMFT)5c + (GW)psc solves for GNon-Locai and 

^ N on — Local 

self-consistently from GNon-Locai = [Gq - 

'^EDMFT — '^Non-Local Gw] ^^„„ r „„„; and Dp 
[Dq^ — ^EDMFT 

^edmft and Uedmft fixed. 

Benchmark - To compare and contrast the approxima- 
tion schemes, we perform a benchmark calculation using 
the model (^) on a 4 x 4, 2D square lattice with peri- 
odic boundary condition. We use, for the free electron 
dispersion, e^: — — (1/2) [cos(/ca;) -I- cos(/cy)]. The half- 
bandwidth is taken as the energy unity. The interac- 
tion is given hy v^:^ U + 2T^[cos(fc2;) -I- cos(A:y)]. We 
study the half-filling case where the strongest correlation 
shows up. We fix the inverse temperature (3 = 8.0 and 



the ratio V/U = 0.25. We vary U from 0.0 to 3.0 for 
the approximation schemes. We benchmark the calcula- 
tion at [/ = 0.5, 1.0, 1.5 by usiiig direct QMC calculation 
via the Hirsch-Fye algorithm |2J|. In all the results we 
are going to present, the major error of the calculation 
comes from the QMC part, in both the exact solution 
and the EDMFT impurity solver. In the latter we solve 
the EDMFT electron-boson irnpurity problem by a hy- 
bridized Monte Carlo method El^lHlll which em- 
ploys an additional continuous auxiliary field. 

At the given V-U ratio, a charge density wave instabil- 
ity at wave vector (tt, tt) is present in the one-shot GW. 
The breakdown is shown in Fig. |3 In the exact solu- 
tion, however, no charge or spin instability is observed 
up to [/ — 1.5. To compare the different schemes on 
the same footing, we restrict ourselves to the paramag- 
netic phase. Depending on the different schemes, this 
may mean to studying a paramagnetic metastable solu- 
tion when the possible instability appears. We find it al- 
lowed for (GW)sc up to the largest U (= 3.0) we studied 
and same for the EDMFT related schemes we described, 
except for (EDMFT+GW)^^ in which the self-consistent 
solution does not exist at C/ = 2.5. 

Our main results are presented in Figs. |3| and 01 
in which the results from the exact QMC and the 
(GW) sc are also plotted. The definitions of the 
plotted quantities are as follows: The Ekmetic and 
Einteraction are the hoppiug and interaction energies 
per site, calculated via the Galitskii-Migdal formula 
■ Eiotai — -^/emetic Ei^iQ^o^ction- The local elec- 
tron self-energy of (EDMFT) 5c is obtained directly 
from the EDMFT solution while for those with spa- 
tial extension, T}°'^{ipn) — (1/iV) I](fc, zp„), with 
TV = 4 X 4. The local Green's function G'°^(ip„) = 
(V^)Efe{[Go(^,«Pn)]"' - T.{k,ip,,)}-\ The quasi- 
particle residue Z = [1— Im I]'°'^(ipo)/po]~^- We measure 
the spatial extension of the electron-self-energy by Zi = 

(l/^)EfeMP(fc,«Pi) -S(fc,ipo)]/(Pi -po)}cosfc,. 

From the overall features of the results, we see: (i) The 
(EDMFT)sc+(GW)psc gives the best results up to the 
largest benchmarked U at 1.5 times the half-bandwidth, 
(ii) The difference between (EDMFT) 5c +(GW)p5c and 
(EDMFT)5c+GW is small quantitatively, since the lo- 
cal self-energies from (EDMFT) 5c is dominant, (iii) 
(EDMFT-l-GW)sc' is not a good scheme in the sense that 
it sees an artificial charge density instability around U ~ 
2.5, where no self-consistent solution is found. Even be- 
fore reaching that regime, the (EDMFT-f GW)5c shows 
significant deviations from the exact solution, (iv) From 
the plottings of Z and Zi in Fig. O we see that (GW)sc 
misses the crossover to localization at large J7, while the 
EDMFT based calculations correctly capture this fea- 
ture. Actually, all the EDMFT related schemes give re- 
sults close to each other at J7 = 3.0 since the spatial 
extensions of the self-energies are no longer important. 
Our results show that, with moderate and strong correla- 
tions, one needs to include higher order contributions be- 
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yond GWA. The schemes, (EDMFT)sc + (GW)psc and 
(EDMFT)5c+GW, offer a reliable solution towards this 
direction. 

Conclusion.- To summarize, we presented a many- 
body scheme which handles the local self-energies non- 
perturbatively via EDMFT and the non-local ones per- 
turbatively via GWA. As an improvement over the lead- 
ing order GWA, the new scheme better captures the 
effects of correlation. We described several implemen- 
tations, which perform self-consistency at different lev- 
els, and benchmarked them by comparing with the ex- 
act solution of a finite-size model system. We found 
that (EDMFT)sc + (GW)psc and (EDMFT)sc+GW 
gave very close results. For the model we studied, 
(EDMFT)sc + (GW)psc gave the best resuh. 

Finally, we should point out that, similar to the self- 



consistent GWA 0, our scheme has the prob- 

lem that the polarization does not have the proper 
asymptotic behavior in the long wave length limit 
[limk^o P {k , uj) oc (fc/cj)^]. Since the schemes combin- 
ing EDMFT with GWA include local vertex corrections 
non-perturbatively, this violation is less severe than that 

in GWA. HilaS 
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